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ABSTRACT 

Several simple experimental acoustic tests of a spray- 
ing system were conducted at the NASA Langley Re- 
search Center. These tests have shown appreciable jet 
noise reduction when an additional cylindrical perme- 
able shell was employed at the nozzle exit. Based on 
these results, additional acoustic tests were conducted 
in the anechoic chamber AK-2 at the Central Aero- 
hydrodynamics Institute (TsAGI, Moscow) in Russia. 
These tests examined the influence of permeable shells 
on the noise from a supersonic jet exhausting from a 
round nozzle designed for exit Mach number, M e =2.0, 
with conical and Screwdriver-shaped centerbodies. The 
results show significant acoustic benefits of permeable 
shell application especially for overexpanded jets by com- 
parison with impermeable shell application. The noise 
reduction in the overal pressure level was obtained up 
to ~5-8’/.. Numerical simulations of a jet flow exhaust- 
ing from a convergent-divergent nozzle designed for exit 
Mach number, M e =2.0, with permeable and imperme- 
able shells were conducted at the NASA LaRC and 
Hampton University. Two numerical codes were used. 
The first is the NASA LaRC CFL3D code for accurate 
calculation of jet mean flow parameters on the basis 
of a full Navier-Stokes solver (NSE). The second is the 
numerical code based on Tam’s method for turbulent 
mixing noise (TMN) calculation. Numerical and exper- 
imental results are in good qualitative agreement. 
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I. INTRODUCTION 

Permeable (perforated) nozzles and other permeable 
devices were studied many years ago for different ap- 
plications, both experimentally and theoretically (see, 
for example, the papers presented at the IAS Meeting 
in Los Angeles, 1954, [1], or the Russian book by G.L. 
Grodzovsky, et al. [2]) . It is well known that compres- 
sion or rarefaction waves in a supersonic flow slowly 
damp when these waves reflect from solid walls or from 
free boundaries (such as density discontinuities). Re- 
search results have shown that it is possible to change 
the supersonic flow structure and flow type by changing 
the reflection quality of the permeable (or porous) wall 
using an appropriate permeability coefficient. In some 
cases, a permeability coefficient can be chosen so that 
the boundary doesn’t reflect incoming disturbances. In 
particular, using such a boundary, we can obtain the 
possibility of controlling the flow velocity (Mach num- 
ber) along the axis of the supersonic divergent nozzle 
portion, of smoothing out any supersonic flow uneven- 
ness, of increasing the permissible size of models tested 
in aerodynamic wind tunnels, etc. 

The use of permeable bodies is a second area of appli- 
cation of this research. Permeability decreases the drag 
coefficient of such bodies because compression waves are 
weakened by reflection from the bow portion of the bod- 
ies with appropriate weakening of the bow shock waves, 
detached or oblique. The simple theory for the solution 
of this problem was proposed by Dr. Kh.A. Rukhmat- 
ulin [3]. The porous and permeable body applications 
are well modeled by the theoretical approach, in partic- 
ular, in boundary condition formulations and by com- 
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parison to the body drag, flow smoothness around the 
body, reduction of noise generated by flow interaction 
with the surface. For example, substantial reduction of 
noise produced by a supersonic jet exhausting from a 
CD nozzle can be obtained by using a porous center- 
body instead of a solid centerbody. This approach was 
proposed in the invention of Dr. L. Maestrello [4], and 
some experimental acoustic test results were presented 
in his paper [5]. 

A third research application connected with perme- 
able surfaces is parachute theory and experiments (see 
[6-8]). In contrast to the previous cases, in the parachute 
theory, different approximations for boundary condition 
formulation use the presence of a substantial normal 
impulse component to the surface which is comparable 
with the transverse component. 

The goal of the research presented is to obtain more 
uniform flow parameter distributions at the nozzle exit 
for reduction of jet noise and simultaneously to increase 
the nozzle thrust. The main mechanism of interaction 
of pressure waves with solid and permeable nozzle wall 
leads to interference of contrary signed waves (compres- 
sion and rarefaction) from the solid wall and lip line por- 
tions. As a result, reflected shock waves interact with 
rarefaction waves and become weaker. Hence, the gas 
flow is more uniform. We have assumed that this effect 
will favour noise reduction generated by the exhaust 
jet. Application of permeable shells or main nozzles, 
in some cases, increases the nozzle thrust. The exper- 
imental and numerical simulation results obtained, in 
general, confirm this assumption. 

II. EXPERIMENTAL APPROACH AND 
ACOUSTIC DATA 

2.1 Experimental tests at the NASA LaRC. 

Several simple experimental acoustic tests for a spraying 
system (Figure 1) were conducted at the NASA Langley 
Research Center. These tests have shown appreciable 
jet noise reduction when additional cylindrical perme- 
able shells were attached at the nozzle exit. Four shells 
were tested, one impermeable and three permeable. All 
cylindrical shells had the same size with the internal 
diameter, d^lin, external, 1 .2in , and length, l=5in; 

the hole diameter was 2^=0.0I5in. 

Two examples of the s hells covering the cylindrical 
pipe of the spraying device are shown in Figures la,b. 
The shell’s location changed during the tests. The main 


tests were conducted using the second shell modification 
shown in Figure lb. The permeable shell has two sets 
of through holes: 8 holes were located uniformly in the 
azimuthal direction along the perimeters of 7 cross sec- 
tions of the cylindrical shell, (i.e. the total hole number 
in each set is 7x8=56). The inclination angle of the hole 
to the cylinder axis a was the same for each set and dif- 
ferent for different sets. The angle changed in the range: 
a — 30°, 45°, 60°, and 90° to the nozzle inlet. At first, 
the impermeable shell was tested. In this case, jet noise 
increased by comparison with noise production with- 
out any shell. Usually, a high frequency tone occured 
during the test with the impermeable shell. Then the 
tests of the permeable shells were conducted and the 
shell location changed by moving it downstream along 
the cylindrical pipe of the spraying device. By this, 
the open hole number increased, which resulted in an 
increase in the mass flow rate through the hole set or 
increase of the permeability coefficient of the shell, e. 
Here s=nx8 xS 0 /S C) where n is number of open hole 
lines (corresponding cross sections with the holes), S 0 
and S c the areas of a round unit hole and of the total 
cylindrical surface of the shell respectively. For some co- 
efficient values, the jet noise was reduced and the high 
frequency tone completely disappeared. 

2.2 Acoustic Tests at the TsAGI, Moscow. 

The reason for these tests was the desire for a more 
accurate examination of this phenomenon. A perme- 
able cylindrical nozzle was made and acoustic tests were 
conducted in the anechoic chamber AK-2 at the Cen- 
tral Aerohydrodynamics Institute (TsAGI, Moscow) in 
Russia. These tests were a continuation and improve- 
ment of some previous acoustic tests using the same 
acoustic measurement apparatus and the same meth- 
ods. The previous test results were presented at the 4th 
AIAA/CEAS Aeroacoustics Conference in Toulouse, 
France, 1998, [9]. We believe that it is expedient to 
repeat some details of the previous approach. 

2.2.1 A Screwdriver centerbody (plug geom- 
etry A Screwdriver centerbody surface belongs to a 
family of shapes formed by rectilinear intervals joining 
corresponding points of two different closed curves in 
space. In this case, a circle as an initial curve and one 
or several crossing rectilinear intervals as an end curve 
are used. Usually, these intervals are symmetrically lo- 
cated relative to a body’s axis of symmetry. 
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One of many possible modifications of such ring-shaped 
nozzles with the Screwdriver and axisymmetric center- 
bodies were made and tested. This design is shown in 
Figure 2 (on the right) and in the large view in Figure 
6. The draft of the meridional plane cross section for 
this design with the axisymmetric (conical) centerbody 
is shown in Figure 4a with the main dimensions of this 
design. 

The 4-petal Screwdriver centerbody (SdCB) shown 
in Figure 6a contains a cylindrical centerbody which 
downstream of the external nozzle exit transfers to a 
Screwdriver shaped portion. There are several geomet- 
ric parameters which define the centerbody: number of 
petals, petal size and centerbody length. The axisym- 
metric centerbody with conical or optimal contour in 
a meridional plane can be taken as a baseline center- 
body foT comparison and for definition of Screwdriver 
centerbody efficiency. In Figure 6a, the 4-petal SdCB is 
based on the conical reference centerbody (CCB) shown 
in Figure 6b. This SdCB surface is formed as follows. 
The 45°-arc of the cylindrical portion and ending ver- 
tical interval are divided by I subintervals [a<, a # '+i... a/ ] 
and [bi,bi+i...bj] respectively, where i=0,l,2,...I. Then, 
the corresponding points a* and are joined by recti- 
linear or curvilinear intervals and these intervals form 
the needed surface. The corresponding cross section ar- 
eas of CCB and SdCB designs are equal. These sections 
are located at the same distance from the nozzle throat 
(or the nozzle exit). 

The manufactured SdCB shown in Figure 6a has curvi- 
linear Screwdriver-shaped surfaces. Namely, longitudi- 
nal rectilinear lines are replaced by curvilinear lines so 
that they are a smooth continuation of corresponding 
straight lines on a cylindrical portion having a horizon- 
tal tangent at the end of a centerbody. Thus, these 
curves should have points of inflection so that they 
may be given by two power functions with conjuga- 
tion at these points of inflection. In the simplest case, 
these functions can be written as follows: Let an ini- 
tial point dj and ending point 6, have Cartesian coor- 
dinates x° Q = 0, y° 0 ~ r 0 cos(&), z° 0 = r 0 sin((/>i), where 
r 0 is the cylindrical portion radius and <f>i is the polar 
coordinate of the initial point. For uniform splitting, 
<f> = (i/I) 7t/ 4, z e = Zg( 1 - i/I). The curve joining these 
points can be described in Cartesian variables as: 


z(0 = z 0 -chZ p ' . 2/(0 = yo+a y t P:> if 0 < 4 < 0 

(1.1) 

*(0 = z e +Mi-O f \ 2/(0 = by( 1-0 73 if Zc < S < 1 

( 1 . 2 ) 

where pi>< 7 i,P 2,?2 are fixed even powers, and £ c is a 
fixed conjunction coordinate of the two power functions 
which can be varied. The coefficients a z } b 2 , a y ,b y are 
defined from the conjunction conditions: equality of the 
function and its first derivatives so that, for example, 
for pi = p 2 — P and qi = q 2 — q these coefficients are: 
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(fo - z t ) 
-p(i - tc )' 


(1.3) 


In accordance with formulae (1.1- 1.3) with p=q=3.0 
and characteristic lengths shown in Figure 4a, nozzles 
with SdCB and CCB were designed and drawn at NASA 
LaRC and then manufactured and tested at TsAGI, 
Moscow in Russia. 


2.2.2 Supersonic CD Nozzle with Permeable 
and Impermeable Shells. The previous acoustic 
tests of the nozzles with the SdCB and CCB have shown 
an absence of any acoustic and aeroperformance bene- 
fits of a SdCB application by comparison with a CCB 
application. Moreover, for some angles 0 = 90°, a noise 
increase was observed ~ 1-2%. The conclusion in the 
report [9] was to continue these tests by moving the 
centerbody into the nozzle and a using permeable noz- 
zle wall. For some reasons, moving the centerbody up- 
stream has been found more difficult than to use con- 
tinuation of the external nozzle cylindrical portion (the 
shell). This shell is mounted on the external main nozzle 
so that the internal shell diameter is equal to the nozzle 
exit diameter and the shell ending axial coordinate, x S) 
is equal to the ending centerbody axial coordinate, x c &. 
The perforated shell has the hole row along the axis 
x containing 23 holes with 3mm diameter. Twenty-six 
such rows are located on the shell cylindrical surface al- 
most uniformly in the azimuthal direction </>. Thus, the 
entire hole number is 23x26=598 which corresponds to 
perforation (permeability) coefficient ~0. 18. All other 
geometric parameters are shown in Figure 4b. 

The nozzle pressure ratio (NPR) is varied in the inter- 
val tt c = 2.5-4.5 with a cold, supersonic jet exhaust. For 
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the nozzle tested with centerbody shown in Figure 4a, 
the design Mach number at the original nozzle exit, i.e. 
at the entrance to the shell, x ey equals ~ A/ e =3.67, and 
nozzle pressure ratio, N PR e = p ( /p e =96.34. These val- 
ues are based on a quasi-one-dimensional theory. There- 
fore, for the nozzle without any shell, the jets are es- 
sentially. overexpanded in the entire interval of total 
pressure, and downstream from the exit intensive shock 
waves are formed. These shock waves reflect repeatedly 
from the centerbody and mixing layer (jet boundary). 
These shock waves are also formed when the solid shell 
is attached to the nozzle. The shock waves reflect from 
the two solid boundaries, centerbody and cylindrical 
wall of the shell. Note that a quasi-one-dimensional 
theory is not accurate enough for this case because of 
the shock waves presence. 

2*2.3 Acoustic Data. The permeable shell can 
weaken shock waves and reduce broadband shock noise 
as well as turbulent mixing noise. For examination of 
this concept, several experimental and numerical tests 
were conducted. Acoustic tests were conducted in the 
anechoic chamber of the Central AeroHydrodynamics 
Institute (TsAGI, Moscow). The interior dimensions of 
the facility within the wedge tips are 9. 6x5. 3x4. 2m high. 
Three nozzle designs with exhaust jets were tested. All 
designs had the axisymmetric external nozzle described 
above and two different internal parts: a) the coaxial 
4-petal Screwdriver shaped centerbody (SdCB), b) the 
coaxial axisymmetric centerbody with the length and 
cross sections equal to the Screwdriver shaped center- 
body (CCB). Three cases were tested: 1) without any 
shell, 2) with perforated shell, and 3) with solid shell. 
All these designs were tested with the two nozzle pres- 
sure ratios: NPR=2.5 and 4.5. In addition, in the case 
lb, an additional test with NPR=3.5 was conducted. 
These cases are enumerated in Table 1 from 1 to 14 for 
more convenient acoustic results illustration. 

The measurement procedure in the anechoic chamber 
and the method of automatic data processing are illus- 
trated in Figure 5. Microphones (model 4136/Bruel& 
Kjaer Co) with cathode followers (model 4633) are po- 
sitioned in the meridional plane on a circular arc with 
the radius i? m =2m with different observation angles 
9 to a positive (downstream) jet axis in the interval 
30° < 9 < 105°. Microphone signals are transmit- 
ted to the magnetic recorder ’’Sony KS-616U” through 


an amplifier (model 2608). Decoding of acoustic pres- 
sure pulsations was conducted by an analyzer (model 
2032/BruelftKjaer Co.) and PC-286 computer which 
provided a narrow band spectrum with a band width 
~ A/— 32Hz. A Pentium computer was used for trans- 
formation of narrow band spectra to 1/3-octave spectra. 

Some experimental acoustic results are illustrated in 
Figures 7-9. In each graphic, two or three curves are 
presented. The curves differ by color and each of them 
are denoted by two numbers divided by a dash (or 
only two numbers for black curves). The first num- 
ber designates the observation angle 9 = 30°, 45°, 60°, 
or 90*\ and the second number designates the case num- 
ber shown in Table 1. Blue lines denote the case without 
any shell, red lines are with perforated shell, and yellow 
lines are with solid shell. 

Acoustic benefits of the perforated shell application 
take place for two observation angles 30° and 45° even 
by comparison with the cases of shell absence. It is re- 
markable because nozzle thrust in the first case is obvi- 
ously more than in the latter cases by additional com- 
pression waves reflected from the shell wall. Numeri- 
cal simulation results confirm this assertion (see below). 
More direct comparison of the two cases, with solid and 
perforated shells, is presented below. Figures 8 and 9 
illustrate substantial acoustic benefits from perforated 
shell application in all observation angles. In Figure 8, 
the narrowband spectral density vs frequency is shown 
for both cases, and in Figure 9, 1/3-octave band spec- 
tral density vs frequency is shown for the same cases. 

Note that in the best cases for the nozzle with con- 
ical centerbody (No. 1-3) and Screwdriver centerbody 
(No. 12-14), the overall acoustic pressure level for the 
nozzle with perforated shell is less than for the nozzle 
with the solid shell up to ~5-8'/*. For the highest noz- 
zle pressure ratio, NPR=4.5, the combination of perfo- 
rated shell and Screwdriver centerbody gives less acous- 
tic benefits (No. 4-6). 

III. NUMERICAL SIMULATION RESULTS 

Numerical simulation of a jet flow exhausting from 
a convergent-divergent nozzle designed with M e =2.0 
and with a permeable shell was conducted at the NASA 
LaRC and Hampton University. Two numerical codes 
were used. The first is the NASA LaRC CFL3D code 
[10] for accurate calculation of jet mean flow parameters 
on the basis of a full Navier-Stokes solver (NSE). The 
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second is the numerical code based on Tam’s method 
[14,15] for turbulent mixing noise (TMN) calculation. 
Numerical and experimental results are in good agree- 
ment. 

3.1 Mean Flow Numerical Simulation. 

Aeroperformance effects and gas dynamic flow char- 
acteristics were analyzed numerically. The nozzle thrust 
calculations were based on a full Navier-Stokes equation 
solver (NSE), and both full and marching Euler codes: 
CFL3D [10], CRAFT [11], and Krayko-Godunov [12]. 
Grid preparation and optimization was conducted us- 
ing GRIDGEN and our own codes. The main results 
were obtained using 2D and 3D versions of the CFL3D 
code which allows the simulation of both inviscid and 
viscous flows. This code is described in detail in the 
CFL3D User’s Manual (Version 5.0) [10]. In accor- 
dance with this manual’s introduction ([10], page 1): 
” ...CFL3D (Version 5.0) is a Reynolds- Averaged thin- 
layer Navier-Stokes flow solver for structured grids... 
CFL3D solves the time- dependent conservation law 
form of the Reynolds- averaged Navier-Stokes equations. 
The spatial discretization involves a semi-discrete finite- 
volume approach. 

Upwind-biasing is used for the convective and pres- 
sure terms, while central differencing is used for the 
shear stress and heat transfer terms. Time advancement 
is implicit with the ability to solve steady or unsteady 
flows. Multigrid and mesh sequencing are available for 
convergence acceleration. Numerous turbulence mod- 
els are provided, including 0-equation, 1-equation, and 
2-equation models. Multi-block topologies are possible 
with the the use 1-1 blocking, patching, overlapping, 
and embedding. CFL3D does not contain any grid gen- 
eration software. Grids must be supplied extraneously.” 

In the paper [9], the coordinate system, grid gener- 
ation, necessary sizes of numerical domains, and mini- 
mal acceptable number of grid points and their distri- 
bution were described in detail. Note here that we were 
used the same as in [9] Menter’s k — w SST turbulence 
model [13]. We will omit other details (refer to the pa- 
per [9] (see p.6 and Figure 1)), and will consider only 
the boundary condition formulation at the permeable 
(perforated) shell wall. In general, that differences the 
presented research from the previous. 


3.1.1 Boundary Condition Formulation for 
Permeable Walls. 

It is well known that compression or rarefaction waves 
in a supersonic flow slowly damp when these waves 
reflect from solid walls or from free boundaries (such 
as density discontinuities). Research results [1,2] have 
shown that it is possible to change supersonic flow struc- 
ture and flow type by changing the reflection quality 
of a permeable (or porous) wall using an appropriate 
permeability coefficient. In some cases, a permeability 
coefficient can be chosen so that the boundary doesn’t 
reflect incoming disturbances. In particular, using such 
a boundary, we can obtain the possibility of controlling 
the flow velocity (Mach number) along the axis of the 
supersonic divergent nozzle portion, of smoothing out 
any supersonic flow unevenness, of increasing permissi- 
ble sizes of tested models in aerodynamic wind tunnels, 
etc. 

The flow at a perforated wall is complicated: through 
some holes, the gas injects into the main flow, and 
through others, a gas suction takes place. These pro- 
cesses become more complicated by viscous effects, in 
particular, by the friction at the solid wall portion. In 
many cases, the detailed flow at the each hole is not im- 
portant, but merely the influence of their presence. The 
solution of the problem for supersonic flow at the per- 
meable wall requires specific boundary conditions. Sev- 
eral approaches are used for such formulations. These 
are enumerated below with a short presentation of the 
final results: 

a) The linearized theory. The assumptions are: 

1. A small difference between pressures at the contrary 
sides of the wall (1 and 2), A p — p 2 —p\ « p, (i=l,2). 

2. A hole size along the wall, X, much larger than the 
hole size in the normal direction to the wall, Y, i.e. 
the velocity component along the wall , u*, is much 
bigger than that in the normal direction v*. 3. The 
wall is plane. The permeable wall is considered as a 
surface with some alternative set of slots, i.e. solid and 
lip shock (density discontinuity) surfaces. Therefore a 
small perturbation theory leads to the simple formulae 
[2]: 

7 = - 1 tf M ' > 1 ( 3 - 1 ) 

Jr = tan Y S )\f l - *7 M ‘ < 1 (3.2) 
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where the main flow at the wall with the corresponding 
Mach number Mi is directed along the wall, x, u , v are 
velocity perturbations in the X and Y directions respec- 
tively, 5 = Sh/($h + S s ) is a perforation coefficient, Sh 
is total hole area, and S s is total solid wall area. These 
formulae provide the boundary conditions which allow 
calculation of the disturbed main flow on both sides of 
the wall. Note that in this approach, any mass flow rate 
through the holes is absent. 


b) The Shock/Rarefaction Wave Method. In 

this approach, a 2D uniform supersonic flow of width 
H, with pressure, p\ , Mach number Mi meets the per- 
meable plane wall portion of length L at the angle, a. 
Then L=H/$ina. The permeability is distributed uni- 
formly along the wall. Dependent upon the pressure 
ratio at the permeable wall, P21 — P2/P1, a shock wave 
(P 2 i > 1) or centered (Prandtl-Meyer) rarefaction wave 
(P 21 < 1) is formed. Thus, instead of at the angle a, the 
flow crosses the wall at some bigger angle 0 = a 4* Ac*. 
The effective permeability coefficient S* is defined as a 
ratio of linear portion of the wall which the flow crosses 
to the wall L*=H/sm0, i.e. S* = L«jL. For pressure 
ratios, P 2 I , close to 1, again simple formulae can be 
written. 

For example, for a cross slot permeable wall: 


7? = A pf 



2(h/b — 1) 

a \J M \ - 1 


(3-3) 


where the effective permeability coefficient 5* — b/h = 
a/(a + Aa) is equal to the ratio between the slot width, 
b, and distance between neighboring slots, h. 


c) The Mutually Penetrating Continua The- 
ory. In accordance with this theory [3,6,8], the perme- 
able surface is considered as a discontinuity over which 
the crossing gas flow loses some normal impulse compo- 
nent, R ny s.t.: 


[p + pu 2 n ] = -Rn, [pu n ll T ] = 0 


(3-4) 


[p?i„] - 0, 


r («’+«?) , « 


+ 


K + 1 p 


*1 = 0 


(3.5) 


where u n , U T are normal and transverse velocity compo- 
nents to the boundary, and k is the specific heat ratio. 
An interphase reaction force R n is defined from special 
experiments or from semiempiric theory. Using a quasi- 
one-dimensional theory inside a hole and the hydraulic 


approximation -[p] = au n + bu 2 , an approximate for- 
mula is presented in the paper [8] for the density ratio 
P 12 = Pl/P 2 : 

P 12 =14- — sin 2 V 5 4- o{M( sin 2 ip) (3.6) 

where £ is a hydraulic loss coefficient depending on the 
permeability coefficient S, local Reynolds number, Re, 
and hole shape. For example, for hole with sharp edges 
and large Reynolds numbers, the £ is: 

£ = - Jj(l -S + 0.707VT^5) 2 
and for holes with smoothed edges: 

i = *(i-s) + (i/.?-i) 2 , *~i 


d) The Porosity Coefficient Definition. In the 

simplest semiempiric theory [1,6], the porosity coeffi- 
cient R p is directly introduced as a ratio between pres- 
sure gradient through the permeable wall, A p and mass 
flow rate through the wall pu n) i.e. 


pJJ 0 u n 

— 7 

A p 


(3.7) 


where U 0 is the imperturbable mean flow velocity along 
the wall, and K p = R p /U 0 is assumed to be a constant 
along the wall. This is for walls made of so called linear 
porous material. For some walls and mean flow with di- 
rection close to normal to the wall, in accordance with 
[6] a quadratic relationship rather than linear is appro- 
priate. The porosity coefficient usually is used as close 
to the permeability coefficient S which is defined above 
as a ratio between hole area and total wall area. 

3.1.2 Numerical Simulation Results. With 
the purpose of understanding and optimization of the 
permeable design for acoustic benefits, several numer- 
ical simulations were conducted at the NASA LaRC 
and Hampton University. Unfortunately, at the present 
time we do not have the capability to calculate noise 
from a 3D jet exhausting from a Screwdriver nozzle. 
Hence, we conducted a numerical simulation for an ax- 
isymmetric CD round nozzle designed for exit Mach 
number, M e = 2.0. The geometric parameters of this 
nozzle are close to the external nozzle tested at the 
TsAGI with axisymmetric and Screwdriver centerbod- 
ies. All the dimensions are shown in Figure 4a. The 
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reference Reynolds number is calculated on the basis of 
critical parameters at the nozzle throat and with the 
characteristic length equal the throat radius, Re * = 
p*r*c„//i„ = 0.128 * 10 5 , where subscripts correspond 
to the parameters at the throat and p is the dynamic 
viscosity coefficient. The cylinder shell (impermeable 
or permeable) is attached to the nozzle exit. The di- 
ameter of the shell is equal to the diameter D e of the 
nozzle exit. The length of the shell is equal to D e . A 
numerical simulation of the jet mean flow and sound 
radiated by the jet have been carried out in this work. 
The mean flow was evaluated on the basis of the NASA 
LaRC CFL3D code. The codes for sound radiation eval- 
uation have been developed by the authors on the basis 
of Tam’s method [14-15]. Three cases have been con- 
sidered: Case I. Overexpanded jet (NPR=6.31). The 
cylinder shell is impermeable. Case II. Overexpanded 
jet (NPR=6.31). The cylinder shell is permeable (coeffi- 
cient of permeability Kp=Cq/Ap=const, Cq is the mass 
flow rate through the cylinder wall per unit area, Ap 
is the difference between pressures on the outside and 
inside walls of the shell, and the value of the nondimen- 
sional coefficient of permeability is chosen to be equal 

0.2). Case III. Underexpanded jet (NPR=9.47). The 
cylinder shell is permeable (Kp=0.2). Some of the re- 
sults of mean flow evaluations are presented in Figure 
10 (Mach contours) and Figure 11 (pressure contours). 
The mean flow evaluation allows us to draw the follow- 
ing conclusions: 

1. The permeable cylinder shell can be used as a 
facility for smoothing of the supersonic jet exhausted 
at the off-design conditions (for both overexpanded jet 
and underexpanded jet). The intensity of barrel-shaped 
shock waves decreases substantially when a permeable 
cylinder shell is used. Barrel-shaped shock waves for 
case I are seen clearly (the upper picture in Figure 10) 
and these shock waves have almost disappeared for the 
cases II and III (the middle and lower pictures in Figure 
10). The cause of the flow smoothness can be explained 
by the fact that pressure along the wall inside the per- 
meable shell changes downstream of the nozzle exit to 
the value of the external pressure (in the ambient air). 
At the same time, for the impermeable shell the jet 
exhausted from the shell is essentially overexpanded. 
(Compare the cases II, III with the case I in Figure 10). 

2. The velocity profile of the flow at the shell exit 
is not uniform. In case II, there is a flux of air into 


the shell through the permeable wall. This flux leads 
to slowing of the flow layers close to the wall. On the 
contrary, in the case III, there is an air flux directed out 
of the shell which causes acceleration of the flow layers 
close to the wall (Figure 11.) 

These jet qualities are important for noise genera- 
tion by the wave emission mechanism as well as by the 
broadband mechanism. In the next paragraph, we will 
discuss the present approach for jet noise calculation. 

3*1.3 Thrust Calculation In accordance with tra- 
ditional thrust definition, introduce P and the corre- 
sponding nondimensional value T as: 

P=f (p e U e 2 + p e )dE- PooH', T = (3.8) 

where subscript indices e, oo, o and * are assigned to 
the nozzle exit cross section, ambient, total and criti- 
cal (in the throat) parameters, £* is the throat area. 
The integration is in some cross section downstream of 
the nozzle exit (or centerbody). The integrand in (3.8) 
is called an impulse function. Such a definition is in- 
troduced for rocket motors, but it does not take into 
account vehicle drag, and assumes the same shape of 
the external and internal vehicle surfaces. Therefore it 
only approximates the real vehicle net thrust. Never- 
theless we will use this definition for an estimation of 
the nozzle shape variation influence on the thrust. 

We can also define the thrust T directly by integra- 
tion of the impulse function at the inlet cross section I Q) 
and the difference between pressure and friction along 
the nozzle wall. The integral of the impulse function at 
the nozzle exit, taking into account the boundary layer, 
allows us to estimate the integral error of the applied 
numerical scheme. Thus the thrust for a single design 
is calculated using the above nondimensional variables 
as: 


r = 5(/ 0 + / 1 )-^.5i, 8 = ^-^—)^ (3.9) 

p 0 7T K + 1 

where 

f r e rl* k 

I 0 — / / p(\ — -M w 2 cjsmot)rdxdrcL<p 

J0 Jr a Jo 2 

(3.10) 


h 



(p + pu 7 )rdrdip 


(3.11) 
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where k is the specific heat ratio, cj is the friction coeffi- 
cient, M w is the local Mach number at the wall, u is ve- 
locity component in the direction of the nozzle-jet axis, 
a is the angle between this axis and the local tangent 
to the wall contour. The correction for calculation of 
thrust losses by the friction effect in the integral (3.10) is 
taken into account only for the cases of numerical simu- 
lation of mean flow in the inviscid approximation based 
on the Euler equations. When the numerical simula- 
tion of the mean flow is based on the full Navier-Stokes 
equations, this term in the integral must be substituted 
by the appropriate shear stress component at the wall, 

t*. 

The calculation of the thrust for all cases of nozzles 
with impermeable and permeable shells based on for- 
mulae (3.8) have produced the following conclusions: 

a) A jet flow exhausting from a convergent-divergent 
nozzle designed for exit Mach number, M e = 2.0 with 
permeable and impermeable shells produces thrust which 
depends on the permeability coefficient K p insignifi- 
cantly for fixed nozzle pressure ratio, NPR, in the in- 
tervals NPR=6.31-9.47 and permeability A" p =0-0.2 con- 
sidered. For example, for the overexpanded jet with 
NPR=6.31, the thrust equals T=0.489, and the devi- 
ation from this value with a variation of permeability 
changes the thrust within the limit of numerical simu- 
lation accuracy for the given grid. But the trust change 
can be significant with increase of nozzle pressure ratio. 
For an underexpanded jet with NPR=9.47 an d K p — 0.1, 
the thrust is equal to 0.456, and for NPR=7.82 (free 
shock conditions) and K p — 0.1, the thrust is equal to 
0.454. At the same time, the influence of the perme- 
ability on the shock wave structure is significant (see 
Figure 10), especially, in the jet. 

b) Similar results were observed for a round noz- 
zle with an axisymmetric centerbody with permeable 
and impermeable shells in the same interval of nozzle 
pressure ratio and permeability. Numerical simulations 
were conducted in the previous work [9] for the nozzle 
shown in Figure 4. a. For this case, the design Mach 
number at the nozzle exit equals A/ e =3.67, i.e. the 
design nozzle pressure ratio is a free shock condition, 
NPR=96.34. Thus, the jets tested are overexpanded 
with the formation of reflected shock waves inside the 
nozzle and shell. The thrust for the nozzle with per- 
meable and impermeable shells varies around the value, 


T=0.420, and also insignificantly depends on the per- 
meability coefficient. 

3,2 Jet Noise Calculation Methods and 
Acoustic Results 

3,2.1 Jet Noise Calculation Methods. 

It is well known ([14]) that turbulent mixing noise 
(TMN) is one of the main components of supersonic 
jet noise in addition to broadband shock noise (BSN) 
and screech tones (ST). In spite of this TMN contribu- 
tion, in a common nozzle, the nearer jet exhaust con- 
ditions come to rated conditions, i.e. when the inter- 
nal jet barrel shock system becomes weaker, the lower 
the noise level produced. The TMN source is turbulent 
small scale pulsations as well as large scale pulsations 
in the jet. With nozzle exit Mach number increase, 
the relative contribution of large scale turbulent pulsa- 
tions to the TMN increases. This phenomenon is due to 
the fact that at nearsonic and supersonic phase speeds, 
sound radiation by these perturbations takes place by 
the mechanism of Mach wave radiation which is very 
effective. Thus, for the case of near perfectly expanded 
high-speed jets, prediction of the jet noise comes down 
to description of the large-scale turbulent evolution and 
evaluation of the sound generated by this turbulence. 
For axisymmetrical jets, a method based on this idea 
has been developed earlier (Tam ft Morris 1980, Tam ft 
Burton 1984, Tam ft Chen 1994). Comparisons of the 
results obtained with the use of this method with the 
experimental data (Troutt ft McLaughlin 1982, Seiner 
at al. 1982) and the results based on direct numerical 
simulation (Mitchell, Lele ft Moin 1997, Mankbadi et al. 
1998) shows favorable agreement. 

The main assumptions underlying the method are the 
following: 

1. Large-scale turbulence is described as a stochastic 
sum of spatially unstable disturbances of the jet (the 
so-called instability waves). Instability waves are the 
disturbances growing in the initial part of the jet ow- 
ing to Kelvin-Helmholtz instability and decreasing at 
downstream parts of the jet where the shear layer be- 
comes thick. The following steps are required to obtain 
the form of the instability waves, a) The mean flow 
in the jet is obtained from experimental data or nu- 
merical simulation of averaged turbulent Navier-Stokes 
equations, b) The non-stationary disturbances (insta- 
bility waves) are obtained on the basis of the solution 
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of linearized Euler equations. The asymptotic method 
of expansion in small parameter and numerical simu- 
lation are used. The small parameter is provided by 
the large disparity in the spatial rate of change of the 
mean flow inside the jet in the radial and axial direc- 
tions. In the first approximation, the spectral equations 
are solved to obtain the most unstable eigen-oscillations 
in the different cross-sections of the jet. In the second 
approximation, the amplitude equation is solved to re- 
late the spectral equation solutions obtained in the first 
approximation for different cross-sections of the jet. 2 . 
Using the method of matched asymptotic expansions, 
the sound radiation generated by separate instability 
waves is obtained. 3* The total turbulent mixing noise 
of the jet is obtained as a stochastic sum of the compo- 
nents generated by separate instability waves. 

At present time, we have developed codes based on 
Tam’s method which can be used for predicting the 
turbulent mixing noise for the case of an axisymmet- 
rical supersonic jet. The next step in this direction 
is a generalization of Tam’s method has be developed 
for the case of jets with cross-section of arbitrary form. 
Such generalization is aimed at jet noise reduction by 
means of nozzle geometry. It is possible that this factor 
may not have significant control over the noise gener- 
ated by fine-grade turbulence (Tam, 1998, [16]). The 
point is that most of the fine-grade turbulence is gener- 
ated in the region near the end of the jet potential core 
where the jet becomes near axisymmetrical independent 
of nozzle geometry. At the same time, the generation 
of the fine-grade turbulence is an essentially non-linear 
process and the maximum turbulence level is indepen- 
dent of the disturbance evolution in the upstream re- 
gions. On the contrary, these reasons do not control 
the large-scale turbulence. This is because the level of 
large-scale turbulence is determined by the increment 
of Kelvin-Helmholtz instability in the initial part of the 
jet. Hence, the noise generated by the large-scale tur- 
bulence should be dependent on the geometry of the 
nozzle. 

3.2,2 Acoustics Calculation Results. Consider 
the influence of mean flow variation caused by the use 
of a permeable shell on the level of the noise gener- 
ated by the jet. As was described above, the supersonic 
jet noise consists of three main components: turbu- 
lent mixing noise (TMN), broadband shock noise (BSN) 


and discrete tones (DT). The components BSN and DT 
are caused by the interaction of non-stationary distur- 
bances and the barrel-shaped shock waves in the jet 
exhausting at off-design conditions. Therefore the level 
of BSN and DT must be decreased together with the 
intensity of the barrel-shock structures when the per- 
meable shell is used. 

The permeable shells also must be useful from the 
standpoint of decreasing of TMN. Indeed, for high-speed 
jets, the main source of the TMN in the direction of 
maximal radiation (30-40 degrees to the jet axis for 
M c = 2) is the large-scale turbulence. In this case, the 
generation mechanism of the TMN is Mach wave ra- 
diation by instability waves. An increase of instabil- 
ity waves is determined by the value of excitement of 
Kelvin-Helmholtz instability. This value has a maxi- 
mum at the initial part of the jet, where the shear layer 
is very thin. It has been noted above that an injection 
of air through the wall of the permeable shell is the 
cause of a slowing down of the layers close to the shell 
wall. In other words, in this case the jet has a relatively 
thick mixing layer, A r m , at the exit of the shell (by 
comparison with the exit radius, r e , namely, A r m /r c 
=0.15). That must diminish the value of the instability 
excitement in the initial part of the jet, and, correspond- 
ingly, must decrease the TMN level. To an even greater 
degree, this conclusion relates to the higher harmonics 
and high-frequency part of the spectrum where the dis- 
turbances have wavelength comparable with the initial 
thickness of the mixing layer. 

This qualitative conclusion is confirmed by the results 
of numerical calculation of the TMN. Some of these re- 
sults are presented in Figure 12. The method permits to 
obtain the results within the accuracy of one arbitrary 
additive constant (in dB scale). This constant is chosen 
to fit the known experimental data on the supersonic 
jet noise (Seiner, McLaughlin, Liu 1982 [17]). The near 
sound field for case II is shown in Figures 12a, b. The 
TMN levels obtained are characteristic of the direction 
of maximum radiation which is approximately at the 
angle of 30 degrees to the jet axis. The comparison of 
the TMN far field levels for the cases I (dashed lines) 
and II (solid lines) at the angle of 30 degrees is pre- 
sented on Figures 12c-e. One can see that an increase 
of the initial value of the mixing layer thickness (from 
0.05 to 0.15) causes a diminishing of the TMN, espe- 
cially for the frequencies St>0.2-0.3 and for harmonics 
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n — 1 and n=2. 

IV. CONCLUSION 

Several experimental acoustic tests were conducted 
in the anechoic chamber AK-2 at the Central Aero- 
hydrodynamics Institute (TsAGI, Moscow) in Russia. 
These tests examined the influence of permeable shells 
on the noise from a supersonic jet exhausting from a 
round nozzle designed for exit Mach number, M e — 2.0, 
with conical and Screwdriver-shaped centerbodies. Sig- 
nificant acoustic benefits of permeable shell applica- 
tions were obtained for overexpanded jets by compar- 
ison with impermeable shell applications. The noise 
reduction in the overal pressure level was obtained up 
to ~5-8'/.. Numerical simulations of a jet flow exhaust- 
ing from a convergent-divergent nozzle designed for exit 
Mach number, Af e =2.0, with permeable and imperme- 
able shells were conducted at the NASA LaRC and 
Hampton University. Two numerical codes were used. 
The first was the NASA LaRC CFL3D code for accurate 
calculation of jet mean flow parameters on the basis of 
a full Navier-Stokes solver (NSE). The second was a nu- 
merical code based on Tam’s method for turbulent mix- 
ing noise (TMN) calculation. The thrust calculations 
for this problem have shown in some cases insignificant 
thrust loss due to permeable shell application and, for 
overexpanded jets, even some thrust augmentation. Nu- 
merical and experimental results are in good qualitative 
agreement. 
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Fig.l Two spraying devices with permeable and impermeable shells tested at the NASA 
LaRC. a) the shell is the long small tickness pipe; b) the shell is the short large tickness 
pipe. 
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Fig.2 The convergent-divergent (CD) conical nozzle (center), Screwdriver and Conical 
centerbodies (right), and solid and perforated shells (left) which were tested in the anechoic 
chamber AK-2 at the TsAGI, Moscow, 



Fig.3 Existing C'D conical nozzle with the Screwdriver centerbody and perforated shell 
mounted in the AK-2 (TsAGI, Moscow). 
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Fig. 4 The picture of the nozzle with the conical centerbody without the shell (a) 
and with the perforated shell (b). The main sizes are in millimeters (mm.). 




Fig. 5 The scheme measurement in the anechoic chamber AK-2 at the 
TsAGI, Moscow, and the scheme of automatic data processing. 


11 





Fig. 6 Axisymmetric convergent-divergent nozzle with the Screwdriver-shaped (a) and ax- 
isymmetric (b) centerbodies. Both centerbodies have the same areas at the cross section 
located on the same distance, x — Xi= const, from the throat, .r = 0. These designs were 
tested in the anechoic chamber at the TsAGI, Moscow. 
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Case No. 

NPR 

Nozzle Configuration 

1 

4.5 

CCB without Shell 

2 

4.5 

CCB with Perforated Shell 

3 

4.5 

CCB with Solid Shell 

4 

4.5 

SdCB without Shell 

5 

4.5 

SdCB with Perforated Shell 

6 

4.5 

SdCB with Solid Shell 

7 

4.5 

CCB without Shell 

8 

3.5 

TTCB without Shell 

9 

2.5 

CCB without Shell 

10 

2.5 

CCB with Perforated Shell 

11 

2.5 

CCB with Solid ShelT^ 

12 

2.5 

SdCB without Shell 

13 

2.5 

SdCB with Perforated ShelT 

14 

2.5 

SdCB with Solid Shell 


Table 1. Acoustic tests at the TsAGI 
Abbreviations: CCB-Conical Centerbody, SdCB-Screwdriver Cenerbody. 




Fig.7a-d. The acoustic power spectral density vs frequency for two observation angles; 
two cases comparison in each figure. 






































NOZZLE WITH IMPERMEABLE AND 
PERMEABLE SHELLS 


MACH CONTOURS 



Fig-10 Mach contours for the supersonic flows inside the CD nozzle designed for the exit 
Mach number, M e = 2.0, and in the exhausting jets for three cases: a) the upper picture- 
impermeable shell (permeability coefficient, K p — 0), overexpanded jet with the nozzle pres- 
sure ratio, NPR=6.31; b) the middle picture-permeable shell with K v — 0.2 and NPR=6.31; 
c) the lower picture-permeable shell with K p =0.2, underexpanded jet with NPR=9.47. 
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NOZZLE WITH IMPERMEABLE AND 
PERMEABLE SHELLS 


PRESSURE CONTOURS 



Fie 11 Pressure contours for the supersonic flows inside the CD nozzle designed for the exit 
Mach number, M e = 2.0, and in the exhausting jets for three cases: a) the upper picture- 
impermeable shell (K p = 0), overexpanded jet with the nozzle pressure ratio, JNPK-b.dl, 
b) the middle picture-permeable shell with K p =0.2 and NPR=6.31; c) the lower picture- 
permeable shell with 0.2, underexpanded jet with NPR=9.47. 
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